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Abstract 

We propose a microscopic model without energy barriers in order to explain 
some generic features observed in structural glasses. The statics can be ex- 
actly solved while the dynamics has been clarified using Monte Carlo calcula- 
tions. Although the model has no thermodynamic transition it captures some 
of the essential features of real glasses, i.e., extremely slow relaxation, time 
dependent hysteresis effects, anomalous increase of the relaxation time and 
aging. This suggests that the effect of entropy barriers can be an important 
ingredient to account for the behavior observed in real glasses. 
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The nature of the glass transition is a longly debated question of much theoretical interest 
It is widely believed that the glass transition is mainly a dynamical process where the 
system can remain trapped in a metastable phase of a finite lifetime depending on the 
rate of the cooling process. If glasses are slowly cooled from the high-temperature region 
then it is possible to reach a crystal phase of very low entropy. But if the system is fastly 
quenched then it reaches a non-equilibrium regime characterized by the existence of very 
slow relaxation phenomena. Usually, the origin of these very slow relaxations is explained 
by the existence of a large number of metastable states separated by energy barriers [[!]]. 
The heights of the energy barriers are widely distributed and the system gets trapped in 
this metastable phase during its time evolution. 

Recently, there have been developments towards a mean-field theory of glasses. In those 
cases one studies systems without quenched disorder with the aid of replica theory . One 
finds the existence of a dynamical transition To where the correlation time diverges. Below 
that temperature the system is always off equilibrium and relaxes towards a dynamical state 
of higher energy. At a lower temperature replica symmetry breaks and a large number of 
states dominate the statics. Below the dynamical transition temperature the system remains 
trapped in a very complex free energy landscape with huge free energy barriers. In mean- 
field theory, the height of these free energy barriers increases exponentially fast with TV 
and metastable states have an infinite lifetime. Generally speaking, free energy barriers 
get contributions from an energetic part and an entropic part. In real structural glasses 
we expect the effects of energy barriers to be substantially different from that in mean-field 
theory because of the existence of nucleation processes 0. Nevertheless, the effect of entropy 
barriers should be not so dependent on the range of the interaction. Then we expect entropy 
barriers to be a relevant mechanism in mean-field as well as in short-range models. 

The purpose of this work is to understand the role of entropy barriers in the behavior 
observed in structural glasses. By entropy barriers we mean the existence of a very small 
number of directions in phase space where the energy decreases. We propose a mean-field 
model with purely entropy barriers and without metastable states. The phase space of 
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this model is very simple and resembles to a golf-hole landscape. It mainly consists of 
flat directions in energy with a very small number of channels where the energy decreases. 
Although the model has no thermodynamic transition it shows a behavior reminiscent of 
real glasses. 

The model. Let us suppose N distinguishable identical particles which can occupy iV 
different states. The Hamiltonian is defined by 

N 

H = -2j2 5 Nr0 (1) 

r=l 

and the energy per site is given by the fraction of occupied states. The N r are the number 
of particles which occupy the state r ( r runs from 1 to N) and they satisfy the constraint 

E N r = N (2) 

r 

The model defined in eq.([l|) can be mapped into a Potts model with a large number of 
states N with Hamiltonian 

N 

H=-J2\m r \ (3) 

r=l 

where m r is the magnetization of state r, m r = N r —1. Looking at the simple Hamiltonian 
of eq.([l]) we observe that there is a trivial ground state with energy per particle ecs — 
—2(1 — -jy) = —2 in the large N limit. In this ground state all the N particles occupy one 
state, its degeneracy being equal to N . The evolution of the system as a function of the 
temperature is as follows. At high temperatures we expect all configurations to have the 
same probability and the energy per particle in this limit is —2(^^-) N = — 2/e. As the 
temperature is decreased the number of occupied states decreases while the occupied states 
increase their occupation numbers N r . Let us suppose we introduce a dynamics for the 
model eq.([l|) (which gives the equilibrium Boltzmann distribution). The time evolution of 
the system as the temperature is decreased is as follows. The rate of variation of the number 
of particles in one state increases due to the particles which reach that state and decreases 
due to the particles which leave that state. The energy decreases when one state is emptied 
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during the dynamical process. Because the total number of particles is conserved, as the 
number of occupied states decreases the time the system needs to empty a further state also 
increases. In this off-equilibrium situation the occupation numbers N r of the occupied states 
perform a random walk and the energy eq.(|I|) decreases very slowly to the static value at 
that temperature. This model has no dynamical phase transition but it shows the onset of 
very slow relaxations in the low temperature region (below T ~ 0.2 close to the maximum 
of the specific heat). We will see that the main characteristics of this model are: strong 
dependence of the energy of the system with the cooling rate, hysteresis effects, anomalous 
increase of the relaxation time and presence of aging. 

Statics of the model. In order to solve the statics of this model we have to compute 
the partition function. We will suppose that we have N particles and each particle % is 
associated with a variable Cj which can take N possible values according to the state the 
particle occupies. 

The partition function is given by, 

1 N 

z = W] T, e M^T, s N r0 ) (4) 

iV • cri r =l 

where the factor N\ in the denominator is a normalization constant in order to make the 
free energy extensive with N. The ocupation numbers iV r satisfy the constraints J2 r N r = N 
and N r = Y^iLi 5 Gi r- Eq.(f|) can be rewritten in term of occupation numbers as 

N i N 

Z= E n N w] exp(2p 1 £5N r0 )5S2N r -N) (5) 

{jV r =0} Hr=l iV »~- r=l r 

We use the integral representation for the delta function, 

5(m) = — / d\e iXm (6) 

Z7T JO 

where m is an integer. Substitution into eq.([5D leads to, 

7 f 2ndX Aif fV exp(i\M + 2p5 M o) ^ m 

z= Jo Y, exvN y- lX+l09 %, wi } ) (7) 

The integral in the previous equation can be readily evaluated by the saddle point method 
in the large N limit. The saddle point A = iz gives the free energy (3f = —Max z A(z) with 



A(z) = -z + log(e w - 1 + e 6 *) (8) 

The saddle point equation is e 2 ^ — 1 = (y — l)e y where y = e z . The solution to this 
equation gives a value y*. The free energy is given by / = —y*//3 and the internal energy 
u = —2e 213 / (y* e y ). We have checked that the first orders in the high-temperature expansion 
of eq.@ for the energy coincide with the previous expression. The energy goes to —2 at 
zero temperature (see fig.l). The specific heat (first derivative of the energy) increases 
approximately like 1/T 2 as the temperature is decreased and shows a maximum at T ~ 0.20. 
At high temperatures the entropy converges to 1 (the number of configurations at infinite 
temperature goes like N N /N\, the dominant term for the entropy being 1). This solution 
is stable but gives a negative entropy below T = 1/(3 ~ 0.345. Because the entropy at zero 
temperature diverges like log(T) (the number of configurations in this limit being N/N\) 
negative entropies are allowed in the model and there is no sign of a thermodynamic phase 
transition. 

Dynamics of the model. We have performed Monte Carlo dynamical calculations in this 
model with the Metropolis algorithm. It is simple enough to allow very fast computations 
for very large values of N. We did two kind of numerical experiments. First we performed 
annealings starting from a random initial configuration. The temperature was slowly de- 
creased and the energy was computed over to MC sweep (a MC sweep corresponds to iV 
trials to change the state of the -randomly chosen- particles). Also we measured the energy 
starting from the ground state configuration (only one state occupied) by slowly increasing 
the temperature. Results are shown in figure 1 for N = 20000 and for the static energy. 
Numerical computations for a larger number of particles (N = 10 5 ) show that finite-size 
effects are negligible. Below T ~ 0.17 we observe a strong dependence of the energy on the 
cooling rate and a slow relaxation of the energy to its equilibrium value. Figure 1 also shows 
the strong dependence of the energy on the temperature change rate during the heating pro- 
cess. The numerical data merges to the static result at a certain temperature. This is also 
the temperature at which the energy departs from the static value in the cooling procedure. 
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The dependence of this merging temperature on the time spent on the cooling rate is an 
estimate of the relaxation time. 

We want now to show that the energy converges to its equilibrium value. We have studied 
the relaxation of the energy at zero temperature starting from an initial random distribution 
of particles. Because there are no metastable states the system can reach the ground state. 
We have measured the time the system takes to reach the ground state at zero temperature 
as a function of the number of particles. We have computed (log{r)) for different values 
of TV ranging from 5 to 20 (the average (...) means average over different random initial 
conditions). We find the typical time very well described by r ~ 0.39 exp(0.67iV). This 
means that the system takes an exponentially large time to reach the ground state. We have 
not succeded in deriving an exact expression for the decay of the energy at zero temperature 
in the infinite N limit. The problem, being highly non trivial, can be approximated taking 
into account the previous result for the exponentially growing time. We argue that the time 
dt the system needs to decrease the fraction of occupied states in a quantity <i(^p) scales 
like 

dt ~ exp(^-)d(^) (9) 

For finite values of iV this expression yields an exponentially large time for reaching the 
ground state. The previous expression means that for a small number of occupied states 
the rate of decrease of the energy is also small (there are less states with more particles per 
state to be emptied). Using u = —2(1 — N oc /N) we get for the decay of the energy 

t = L duexp ^2T^ (10) 

where Uq is the initial energy at time zero. We have measured the decay of the energy 
as a function of time. While this expression is only approximate it shows a remarkable 
agreement with the numerical data especially in the large time region over several decades 
of time. The relaxation of the energy, far from being of a logarithmic or algebraic type, is 
extremely slow as eq.(|II]) shows. The fact that the energy converges to the static result at 
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zero temperature suggests that a dynamical transition at a finite temperature is absent. In 
what follows we will check this point by computing the relaxation time. 

The relaxation time and aging. To fully characterize the dynamics of this model we have 
computed the relaxation time. We can define two types of correlation functions, one for the 
(jj variables, the other one for the energy state variables. In the regime of low temperatures 
the system performs a random walk changing particles from one state to another and the 
appropiate correlation function is given by the state to state energy function, 

r l, ^ - gkM!> - "(<'))M0 - "(0) ,,,, 

Ce(M) s£.,Mt) -«(«))> ( ' 

where e r (t) = 5jv r (t)o an d u(t) is the mean energy per site at time t. We have normalized 
it in order to have C e (t, t) = 1. For times t larger than the correlation time t(T) the C e (t, t') 
should depend only on the time difference t' — t (t' < t) and decay exponentially with time 
(C e (t,f) ~ exp(—(t' — t)/r)). For t « r(T) the system is off-equilibrium, aging efects are 
present and time translation invariance is broken. We have measured the relaxation time as 
a function of the temperature in the low temperature region. Calculations for this model are 
fast enough to allow to compute several orders of magnitude in the relaxation time. Results 
are shown in figure 2. There is no sign of algebraic divergence at any finite temperature. 
While a pure Arrhenius divergence r ~ exp(A/T) does not fit enough well (data present a 
systematic curvature in figure 2) we find that a Vogel-Fulcher law r ~ exp(A/(T — T )) 0] 
describes extremely well the increase of the relaxation time. The value of T ~ 0.02 is stable 
to including more points in the fit and is definitely better than the Arrhenius one. We do not 
attach special physical meaning to the value of To but the indication of an anomaly in the 
divergence of the relaxation time. This is different to the case of models where metastability 
is present where T can be identified as a thermodynamic transition temperature ]|]. 

The dynamical transition in this model takes place at T = and we expect for t finite 
the correlation function of eq.(|ll|) to display aging. Introducing the waiting time t w and 
redefining the times t w = t, t' = t + t w in eq.flTTD we find that C e (t w ,t w + t) is pretty well 
described by the scaling law 



C e (t w ,t w + t) = f(-) (12) 

Similarly to other mean-field models ||[7|, it is very plausible that this scaling behavior 
is exact in the large t w limit at least at zero temperature. Results are shown in figure 3. 
Data collapse in a single curve, the scaling function / of eq.(0) scales like f(x) ~ 0.78 x~ ' 70 
for large values of x. If we define uea = limt u ,-+<» C e (t w , 2t w ) we find that this value jumps 
discontinuously to a finite value ~ 0.58 at T = being zero at any finite T in agreement 
with the absence of any finite T dynamical transition. 

We can summarize now our results. We have introduced a simple model without energy 
barriers (i.e. without metastable states) and without disorder. We have exactly solved the 
statics while the dynamics has been clarified by numerical computations of the model. Re- 
laxation is extremely slow at low temperatures due to the presence of high entropy barriers, 
i.e. existence of small number directions in phase space where energy decreases. While the 
system has no thermodynamic transition it displays the main features of real glasses, i.e. 
extremely slow relaxation, time dependent hysteresis effects, anomalous increase of the re- 
laxation time and aging. This suggests that the presence of activated energy barriers are not 
the only possible ingredient in order to find a Vogel-Fulcher behavior (and eventually Arrhe- 
nius behavior) for some transport quantities. The fact that the relaxation time can be nicely 
fitted to a Vogel-Fulcher law in this model is an indication that the existence of anomalies 
in the relaxation time can be affected by entropy barrier effects. The behavior observed 
in this model suggests that it is not easy to disentangle the effects of energy barriers form 
the effects of entropy barriers, at least experimentally. Real glasses do have energy barriers 
and it is clear that this model cannot explain all the experimentally observed features (for 
instance, existence of a cristallization transition). Interestingly enough, it presents the main 
dynamical features observed in structural glasses. The presented model is of mean-field type 
because a spatial arrangement of the states is absent and it is possible to include spatial 
correlations. Also one can consider generalizations of this model, for instance, defining the 
energy as given by the number of occupied states above a certain finite level k. In this last 
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case we expect to find similar results to those presented here. 
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Figure Captions 

Fig. 1 Energy as a function of temperature for two different processes (cooling and heat- 
ing) and three different cooling rates t (defined as the number of MC sweeps per 
temperature step, this being 0.005). The continuous line is the static energy. 

Fig. 2 Relaxation time as a function of l/T. The continuous line is the Vogel-Fulcher law 
r = Aexp{B/(T - T ) with parameters A = 0.89,5 = 0.51, T = 0.02. The dotted 
line is the fitted Arrhenius law r = Aexp(B /T) with A = 0.63, B = 0.66. Error bars 
for the relaxation time are of the size of the crosses. 

Fig. 3 C e (t w ,t w + t) for iV = 100000 as a function of t/t w for different t w values at zero 
temperature. 



11 



10 6 
10 5 
10 4 
1000 
100 
10 



1 





0,8 




tw+t) 


0,6 


A 

O 






+ 

X 

□ 
♦ 


Ce( 


0,4 




0,2 


o 

ffl 



tw 

30 

100 

300 

1000 

3000 

10000 

30000 

100000 



o 1 

1 0" 4 1 0" 3 1 0" 2 



10° 10 1 10 2 



